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The Fermi Gamma-ray Space Telescope, as an all-sky survey and monitoring mission, is producing daily/weekly 
sampled gamma-ray light curves for dozens of blazars and other high-energy sources. Highlights on MeV-GeV 
gamma-ray variability properties of these sources are reported together with a few remarks about some multi- 
waveband observing campaigns led by Fermi and targeted to known or newly discovered gamma-ray blazars. 



1. Gamma-ray MeV-GeV band variability 
of Fermi blazars 

The Large Area Telescope (LAT), on board the 
Fermi Gamma-ray Space Telescope Atwood et al. 
[2009], is a pair-conversion 7-ray telescope, sensitive to 
photon energies from about 20 MeV up to > 300 GeV. 
The LAT consists of a tracker (two sections, front and 
back), a calorimeter and an anti-coincidence system to 
reject the charged-particle background. Fermi LAT, 
working in all-sky survey mode, is an optimal hunter 
of high-energy flares, transients and new gamma-ray 
sources, and is an unprecedented monitor of the vari- 
able 7-ray sky, thanks to the large peak effective area, 
wide field of view (~ 2.4 sr), improved angular reso- 
lution and sensitivity. 

The entire 7-ray is observed every 2 orbits 3 
hours) representing a continuous monitor for variabil- 
ity and transients allowing the collection of regular 
daily /weekly-sampled light curves for dozens of GeV 
sources. Multiwavelength (MW) observing campaigns 
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Figure 1: Multi- waveband (a.k.a. multiwavelength, MW, 
multifrequency) variability observed from radio-band to 
TeV 7-ray s. Yellow shaded polygon roughly represents 
the timescales range probed for variability at different 
energy bands (i.e., measures in the power density spectra 
versus the spectral energy distribution space, PDS-SED 
plane). The Fermi LAT instrument is exploring a large 
portion of this plane. Adapted from Wagner [2001]. 
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are therefore limited only by the source brightness and 
by the ability to coordinate other telescopes (Fig. 1). 
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Figure 2: Top panel: daily-bin 7-ray flux light curve of 
blazar 3C 454.3 (100 MeV - 200 GeV band, red/gray 
points, main panel) between MJD 55070 and 55307 (2009 
Aug. 27 2010 Apr. 21). Dashed black lines mark period 
over which the SF, PDS and Wavelet analysis are 
conducted (using improved resolution with 3 hour bins). 
The light curve of the previous 2008 JulyAug. flare, 
shifted by 511 days, is also shown (violet/dark-grey tiny 
points). Insets show blow-ups of the two periods (A and 
B on the plot) when the largest relative flux increases 
took place (red/gray filled, blue/dark open, 
green/light-gray filled data points corresponding to daily, 
6- hour, 3- hour bin fluxes, respectively). Bottom left 
panel: SF of the 3h-bin flux light curve for the period 
MJD 55140-55259 (2009 Nov. 5 - 2010 Mar. 4, black 
dashed lines in the top panel) and corresponding PDS 
(inset). Bottom right panel: 2D contour plot of the 
continuous Morlet wavelet transform power density of the 
same light curve (thick black contours: 90% confidence 
levels of true signal features against white-noise and 
red/flickering noise background; cross-hatched region: 
cone of influence dominated by spurious power of 
time- frequency edges). From Ackermann et al. [2010]. 
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Figure 3: Main panel: 3 1.5- month (945-day) light curve of the integral 7-ray flux above 200 MeV in weekly time bins 
of blazar PKS 1830-211, from 2008 August 04, to 2011 March 7 (MJD 54682.65 to 55627.65). This is an example of the 
Fermi LAT capabilities in high-energy temporal variability monitoring. Inset panels: 12-hour bin light curves detailing 
the period around the mild flare of October 2009 (A interval), detailing the period around the large outburst of 
October 2010 (B interval) and the secondary double flare of December 2010 and January 2011 (C interval). 




Figure 4: Results from 7-ray light curves extracted in fixed 1-month bins over the first 2 years of Fermi LAT all-sky 
survey for the 886 blazars/AGNs of the clean source sample of the second Fermi LAT Source Catalog, 2FGL. Left 
panel: variability index versus isotropic 7-ray luminosity (red/grey: FSRQs; blue/black: BL Lac Objects). Dashed line 
represents the 99% confidence level for a source to be variable. Central panel: distribution of the temporal PDS 
power-law indexes (a = /3 + 1) for the FSRQs (red/continuous line) and BL Lac objects (blue/dashed line) of the 
sample evaluated in time domain using first order structure function (SF) analysis (blind power-law index /S estimation 
using a maximum lag of 2/3 of the total light curve range). Cumulative power density spectra (PDS) for bright FSRQs 
(red/top line) and BL Lac objects sub-samples (blue/bottom line) showing similar slopes (inset). From Ackermann et 
al. [2011]. Right panel: scatter plot of the weekly-binned variability auto-correlation (DACF) crossing times in the rest 
frame of the source (corrected for z and beaming) versus the total apparent isotropic 7-ray luminosity {E > 300 MeV) 
in the co-moving frame for 15 bright LAT blazars that are also part of the MOJAVE program. 3C 454.3 is out of the 
plot range (with At^esf^ = 254.3) weeks and log^Q(L£;) =48.1. Adapted from Abdo et al. [2010a], Ackermann et al. 
[2011]. 



Irregular and aperiodic variability is found in 
blazars at all the timescales and at all the ener- 
gies (MW variability). EGRET already showed that 
blazars have a high-energy component in their spec- 
tral energy distributions (SED) and are the largest 
class of variable 7-ray sources, although it was Hm- 
ited by statistics, while Fermi is observing low 7-ray 



brightness states too (example: BL Lac in Fig. 6). 
Of the studied Fermi LAT blazars, 2/3 are variable 
and high states are less than 1/4 of the total light 
curve range [Abdo et al. 2010a]. Flat Spectrum Ra- 
dio Quasars (FSRQs) and low energy peaked BL Lac 
objects (BL Lacs) show the largest relative variance 
(Fig. 4), while high-energy peaked BL Lacs display 
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Figure 5: Radio, mm, optical (flux, polarization), UV, X-ray, 7-ray fluxes light curves obtained by Fermi coordinated 
MW campaigns (PKS 1502+106, left, 3C 279, right). Adapted from Abdo et al. [2010b], Abdo et al. [2010c]. 



a lower variable but persistent emission. Sources like 
PKS 1510-08, PKS 1502+106 (Fig. 5 and 6), 3C 454.3 
(Fig. 2), 3C 279 (Fig. 5), PKS 1830-211 (Fig. 3), 4C 
21.35, 4C 38.41 (all FSRQs) and AO 0235+164, 3C 
66 A (both BL Lacs) are among the brightest, most 
variable, and isotropically luminous blazars seen by 
Fermi . For bright flares even intra-day light curves 
have been extracted as can be seen, for example, in 
Fig. 2 and 3. 

Discrete autocorrelation function (DACF) and 
structure function (SF, Fig. 4 and 2) analysis showed 
different patterns, autocorrelation times and power- 
law PDS indices (1//" fluctuations, where / = 1/t) 
implying different variability modes for each source 
(more flicker, a 1 or more Brownian, a > 2, dom- 
inated). 3C 454.3 is a fully Brownian 7-ray source 
while other powerful blazars have values half-way be- 
tween the two modes. Average power spectral density 
analysis (PDS) in frequency, /, domain over the blazar 
subclasses points out a slopes from 1.3 to 1.6 [Abdo 
et al. 2010a, Ackermann et al. 2011]. No evidence for 
persistent characteristic 7-ray timescale(s) is found. 
Flare profiles are mostly symmetric (in part because 
of superimposing flares and large bin smoothing). The 



fractional variability during outburst appears similar 
to its longer term mean in the few objects studied in 
detail. 



2. Multi-waveband blazar variability 
through observing campaign led by 
Fermi 

Fermi-driven MW observing campaigns are shed- 
ding light on the PDS-SED plane (i.e. timescale- 
energy parameter space). Broad-band MW studies 
are addressed mostly to cross-correlation and time- 
lag analysis, to time-resolved SED modeling, to the 
search for orphan flares and spectral hysteresis, to 
the analysis of 7-ray vs synchrotron amplitude ratios 
and emission peaks and to the study of the radio-7- 
ray connection and source populations. Simultaneous 
MW observations are crucial also for the identification 
of newly discovered 7-ray sources. 

In order to define and better constrain physical pa- 
rameters, processes and emission components, to clar- 
ify the role of the central engine, jets and their in- 
terplay, the jet composition and structure in AGNs 
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Figure 6: Multi-epocli spectral energy distributions of 
blazars (PKS 1502+106, BL Lac, 3C 66A) with 
multifrequency data obtained with coordinated and 
simultaneous MeV-GeV gamma-ray (by Fermi LAT) and 
radio, mm, optical UV, X-ray and TeV (bottom panel 
only) observations. Adapted from Abdo et al. [2010b], 
Abdo et al. [2011b,c]. 



and blazars is necessary to collect more different and 
longer sequences of MW observations. Some clues are 
already emerging. 

1) The knowledge of redshifts is crucial but ~ 50% 
of BL Lacs have still unknown z. 2) Simple single- 
zone synchrotron self Compton (SSC) descriptions are 
vanishing. 3) Internal shock scenario, eventually with 
composite particle energy distributions, works well. 
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4) Cross-correlation analysis with optical polarization 
provides important clues on jet physics. 5) The loca- 
tion of emission site can be both inside and outside the 
broad line region (BLR). 6) Magnetic fields are com- 
\ plex but can be highly ordered and jet-aligned during 
'„ 7-ray flares. 7) In some cases bright 7-ray flares seem 
i to occur after ejections of superluminal radio knots. 
^ Two examples of composite MW light curves obtained 
Z from Fermi blazar campaigns (PKS 1502+106 and 3C 
279) are reported in Fig. 5, Abdo et al. [2010b], Abdo 
et al. [2010c]. 

In conclusion Fermi LAT, in the frame of blazars 
and other AGNs science, is demonstrating very good 
capabilities in the field of 7-ray variability analysis 
and radio-gamma-ray connection and is showing an 
optimal synergy with the Swift mission. 



Acknowledgments 

The Fermi LAT Collaboration acknowledges support from a 
number of agencies and institutes for both development and the 
operation of the LAT as well as scientific data analysis. These 
include NASA and DOE in the United States, CEA/Irfu and 
IN2P3/CNRS in France, ASI and INFN in Italy, MEXT, KEK, 
and JAXA in Japan, and the K. A. Wallenberg Foundation, 
the Swedish Research Council and the National Space Board in 
Sweden. Additional support from INAF in Italy and CNES in 
France for science analysis during the operations phase is also 
gratefully acknowledged. 



References 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010a, 

ApJ, 722, 520 
Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010b, 

ApJ, 710, 810 
Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010c, 

Nature, 463, 919 
Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011a, 

ApJ Lett, 733, L26 
Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011b, 

ApJ, 730, 101 
Abdo, A. A., Ackermann, M., Ajello, M., et al. 2011c, 

ApJ, 726, 43 

Ackermann, M., Ajello, M., Baldini, L., et al. 2010, 

ApJ, 721, 1383 
Ackermann, M., Ajello, M., Allafort, A. et al., 2011, 

ApJ, accepted, (arXiv: 1108. 1420) 
Atwood, W. B., Abdo, A. A., Ackermann, M., et al. 

2009, ApJ, 697, 1071 
Wagner, S. J. 2001, ASP Conference Series, 227, 112 



